A microscopic view of accelerated dynamics in deformed polymer glasses 
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A molecular level analysis of segmental trajectories obtained from molecular dynamics simula- 
tions is used to obtain the full relaxation time spectrum in aging polymer glasses subject to three 
different deformation protocols. As in experiments, dynamics can be accelerated by several orders of 
magnitude, and a narrowing of the distribution of relaxation times during creep is directly observed. 
Additionally, the acceleration factor describing the transformation of the relaxation time distribu- 
tions is computed and found to obey a universal dependence on the global strain, independent of 
age and deformation protocol. 
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When amorphous polymers are being deformed, the 
slow glassy dynamics resulting from broad distributions 
of relaxation times becomes accelerated and permits plas- 
tic flow 0, Q • Recent experimental advances have made 
it possible to directly measure segmental relaxation dy- 
namics in poly(methyl methacrylate) through an optical 
photobleaching technique that probes the reorientation 
of small dye molecules 0, 0|- During deformation un- 
der applied load, molecular mobility increases by up to 
three orders of magnitude with a concomitant narrowing 
of the relaxation time spectrum Hi]. Additionally, shifts 
in mobility relative to the undeformed glass strongly 
correlate with the instantaneous strain rate. In order 
to explain plasticity in amorphous solids, the literature 
frequently invokes Eyring's phenomenological picture of 
stress-biased activated processes, in which the viscosity 
and hence relaxation times 7y decrease with applied stress 
a as r r ~ a I sinh [aV/kBT] Q. It has been demonstrated 
multiple times, however, that this simple model fails to 
predict the yield behavior of polymer glasses consistently 
over a range of stresses, strain rates and temperatures [6j. 
For instance, relaxation times change even when the ap- 
plied stress is constant 

Enhanced mobility in deformed polymer glasses has 
also been observed in both atomistic and coarse grained 
computer simulations. Chain dynamics was probed by 
measuring the monomer mean-squared displacement [7|, 
the dihedral transition rate between trans/gauche states 
Q , and the decay rate of covalent bond autocorrelations 
and intermediate scattering functions 0, [l(| ■ These very 
different observables all report increased local mobility 
on the intra- and interchain level both in the sub-yield 
Q and post-yield Q regime and confirm the experimen- 
tally observed strong correlation between relaxation time 
and strain rate [ij [lfj. A limitation of these studies is, 
however, that they only report spatially averaged quan- 
tities and do not reveal the full spectrum of relaxation 
times. Furthermore, since even undeformed glasses are 
not in equilibrium, it is important to consider the inter- 
action between deformation and physical aging. Aging 



increases the stiffness of glassy polymers by increasing 
structural relaxation times; however, deformation alters 
the intrinsic aging dynamics, and plastic flow can erase 
the history of the glass, an effect known as "mechanical 
rejuvenation" [ll|, |T2| ■ 

In the present Letter, we use molecular simulations as a 
computational microscope to obtain detailed insight into 
the dynamics of polymer glasses under deformation, and 
its consequences for physical aging. From an analysis 
of individual particle trajectories, we identify relaxation 
events and hence obtain the full probability distribution 
of relaxation times. We consider three different defor- 
mation protocols: a step stress (creep), a constant strain 
rate deformation, and a step strain (stress relaxation). 
The relaxation time distributions are changed in complex 
ways which depend on the wait time and the deformation 
protocol, however we find that these changes can all be 
accounted for using a universal acceleration factor which 
depends only on the global strain. 

Our study employs a well-known coarse-grained bead- 
spring polymer system, which has been used extensively 
to model glassy dynamics flij and polymer deforma- 
tion [lij]. The beads are bonded via a finite extensi- 
ble non-linear elastic (FENE) spring and interact via a 
truncated 6-12 Lennard Jones (LJ) potential Vi,j(r) — 
4uo [(a/r) 12 — (a/r) 6 ] for r < r c — 1.5, where uq and a 
set the reference energy and length scale. Results will 
be given in LJ units, where the reference time scale is 
tlj = y ma 2 Juq. Individual samples consist of 20,000 
monomers forming short chains of length 10 in a cubic 
box subject to periodic boundary conditions. After equi- 
librating at a melt temperature T = 1.2, the glass is 
prepared through a rapid quench at constant volume to 
a glassy temperature of T = 0.25, where the hydrostatic 
pressure is close to zero. The glass is then aged at zero 
pressure for a waiting time t w before a uniaxial deforma- 
tion is applied in the form of (i) a step stress of magnitude 
cr, (ii) a constant strain rate e, and (Hi) a step strain of 
magnitude e. In all three protocols the stress components 
perpendicular to the deformation axis are kept at zero, 
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FIG. 1: (Color online) Distributions of first hop times t\ (a)-(c) and persistence times r (cl)-(f) in an aging polymer glass 
subject to deformation after undergoing physical aging for three waiting times t w = 750 (black), 7500 (red/gray), and 75000 
(green/light gray). Panels (a) and (d) correspond to a step stress with amplitude a = 0.4 (dashed) and 0.5 (dotted). Solid dots 
show p(r) for t < r a (see text) and t w = 75000. Panels (b) and (e) refer to deformation at imposed strain rate e = 8.9 x 10~ 7 
(dashed) and 8.9 x 10~ 6 (dotted), while panels (c) and (f) show the effect of a step strain with amplitude e = 0.02 (dashed) and 
0.04 (dotted) For comparison, the corresponding distributions in an undeformed glass that is only undergoing physical aging 
are also shown in every panel (solid lines). Straight lines indicate power law with the given slopes. Insets show the mechanical 
response to (a) a = 0.4 (solid) and 0.5 (dashed), (b) e = 8.9 x 10 -7 , and (c) e = 0.04 for the three wait times. 



while stress/strain/strain rate are maintained at the in- 
dicated amplitudes. 

In amorphous solids, the motion of individual particles 
is not smooth, but highly intermittent. Atoms spend long 
times in caged environments before undergoing a rapid 
structural relaxation to a new position. These relaxations 
correspond to hops on a segmental trajectory. We have 
developed an algorithm that reliably identifies such hops 
through peaks in the standard deviation of the particle 
position (see ref. for details) . The analysis yields the 
full temporal sequence of relaxation events for a subset of 
5000 polymer beads chosen randomly from the full sim- 
ulation volume, from which relaxation time distributions 
are readily constructed. 

In Fig. [T] we show the observed distributions of hop 
times for the three deformation protocols defined above 
and the undeformed glass. We separately present the 
distribution of first hop times p(t\) since the onset of de- 
formation (panels (a)-(c)) and the distribution of times 
in between all subsequent hops p(r) (panels (d)-(f)). r 
measures the lifetime of an individual caged configura- 
tion, which we call the persistence time. The mechanical 
response of each protocol is shown as an inset in pan- 



els (a)-(c). In ref. [lfjj, we showed that for simple aging 
with no deformation, the first hop times increase with 
a power law in the wait time, t\ ~ tfo, where fj, is the 
aging exponent. The persistence time distribution, how- 
ever, is age independent and takes the form of a power 
law p(r) ~ -r -1 - 2 . Since the mean hop time of such a dis- 
tribution is infinite, the dynamics is non-stationary and 
ages. Studying the changes to these distributions un- 
der deformation provides rich insight into the origin of 
deformation-accelerated dynamics. 

First, we consider the effects of a constant stress de- 
formation. The inset of Fig. Ufa) shows that the strain 
increases with a characteristic timescale r a , which grows 
with wait time as tfa. In the sub-yield regime, the re- 
sponse remains shifted in time for different wait times, 
whereas above yield, the glass eventually begins to flow 
and all aging is erased (mechanical rejuvenation). The 
first hop time distribution p{t\) (Fig.QJa)) is significantly 
narrowed for increasing stress. Not only is the probability 
to observe short first hop times increased relative to the 
undeformed sample, but the power law tail decays more 
rapidly. In this way, our analysis permits direct obser- 
vation of the narrowing of the relaxation time spectrum 
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FIG. 2: (Color online) Strain rate versus hop rate for a — 0.5 
(o), 0.4 (□), 0.3 (A), 0.2 (0) and t w = 750 (black), 7500 
(red/gray), and 75000 (green/light gray). 

during creep deformation, which previously was inferred 
only indirectly from stretched exponential fits to correla- 
tion functions [4(. 

Fig. [TJd) shows the persistence time distributions, 
which for a constant stress deformation depend explicitly 
on both wait time and the time since the deformation be- 
gan (measurement time t) . Lines in Fig. QJd) show p(r) 
for cages that come into being after the alpha relaxation 
time (t > t q ), when the strain is increasing slowly. Sim- 
ilar to the first hop time distribution, p(r) is narrowed 
and decreases faster for large r. This effect is even more 
pronounced when p(r) is calculated from particles whose 
cages were created at t < r a (dots) . The narrowing of the 
persistence time distributions provides an interesting ex- 
planation for mechanical rejuvenation. Our results show 
that at a = 0.5, where significant plastic flow occurs, the 
tails of p(r) decay with an exponent of ~ —2.2. For such 
a distribution, the mean persistence time becomes finite 
and a steady state is reached which does not depend on 
the history of the sample. 

Understanding the dynamics during creep has been the 
focus of several recent investigations [3, 0, E3] • In these 
works, the mobility is defined as the inverse of the average 
relaxation time, computed through repeated measure- 
ments of an autocorrelation function during the course of 
the experiment. Results show that the mobility increases 
almost linearly with the strain rate, independent of stress 
and temperature. The mobility as defined in these ex- 
periments is conceptually similar to the total hop rate, 
accessible through our technique. In Fig. ® the strain 
rate is plotted versus the hop rate during the step stress 
experiment. Results for all stresses and wait times fall 
on the same, universal curve. The curve is not exactly 
linear due to the presence of a finite hop rate even in the 
undeformed, aging samples, and an increase of the mean 
(non-afiine) displacement per hop with strain rate. 
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FIG. 3: (Color online) (a) Cumulative probability distribution 
of the first relaxation event for the undeformed glass (solid 
line) and under deformation with a step stress a = 0.4 (o), 
a step strain e = 0.01 (A), and a constant strain rate e = 
8.5 x 10" 7 (□) for t w = 75000. (b) Acceleration ratio t u /t d 
as a function of global strain e for three different deformation 
protocols. Stress step: a = 0.3 (o), 0.4 (□), 0.5 (A); constant 
strain rate: e = 8.9 x 10~ 6 (0), 8.9 x 10" 7 (V); strain step 
(<). For each: t w = 75000 (black), 22500 (red). 



A constant strain rate deformation produces a very dif- 
ferent effect on the hop times than the creep experiment. 
Now the system is undergoing plastic flow at an exter- 
nally imposed rate, and this timescale is reflected in the 
relaxation time distributions. One can see in Fig. [lib) 
that the tails of the p(ii)-distributions steepen consider- 
ably at times of order 0.1/e. At short times, however, the 
distributions are unmodified relative to the unperturbed 
glass. Similarly, the persistence time distributions p(r) 
(Fig. [He)) are not modified at short times, but decrease 
rapidly over the same timescale as p{t\). p(r) remains 
independent of the wait time and the measurement time. 
Aging is also erased in this protocol, but now due to trun- 
cation of the power law distribution rather than a change 
in exponent as was observed for constant stress. These 
results are completely consistent with the effect of aging 
on the macroscopic mechanical behavior (inset): stress- 
strain curves display an initial overshoot stress that in- 
creases with age, followed by an age independent flow or 
hardening regime. 

Finally, we analyze the hop statistics for the step strain 
experiment. Fig. (He) shows that increasing the am- 
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plitude of the strain step primarily modifies the short 
time part of p(t±). Alternatively, the persistence times 
(Fig-Elf)) are completely unchanged relative to the unde- 
formed glass. The effects of a strain step on the first hop 
time can more easily be appreciated through the cumu- 
lative distributions P(t x ) = f* 1 p(*'i)^'i- 1 - P(t) is the 
probability that a particle has not yet hopped at time t, 
and its decay closely resembles the autocorrelation func- 
tions more typically used to measure particle relaxations 
such as the intermediate scattering factor 17|. Figure 



|3fa) shows that the cumulative distribution in the step 
strain protocol is simply shifted in time by a constant fac- 
tor with respect to the undeformed glass: all relaxation 
times are rescaled identically. 

The cumulatives for the step stress and the constant 
strain rate protocol are also shown in Figure J3ja). Un- 
like the simple time shift seen in the step strain exper- 
iment, undeformed and deformed curves continue to di- 
verge with time. To make the transformation of the re- 
laxation time distributions more quantitative, we define 
an acceleration ratio as the ratio of times when the unde- 
formed and deformed cumulative distributions are equal: 
ti/tf, where t" and tf are defined by the relationship 
P(ti) = P(tf). The acceleration ratio captures changes 
to the full distribution of relaxation times and is there- 
fore distinct from the average mobility defined in earlier 
works. Remarkably, the variation in acceleration ratio 
for the different deformation protocols can be reduced 
to the effect of their global strain. A parametric plot 
of the ratios t\jt\ versus the total strain experienced 
by the sample is shown in Fig. EJb) for all three pro- 
tocols and two different waiting times t w [lij . Despite 
the dramatically different mechanical response, all data 
sets collapse reasonably well on a common curve of the 
form tf/tf ~ exp(/(e)), where / is a universal function of 
strain. We were unable to achieve a similar collapse us- 
ing any other deformation variable. For sub-yield strains, 
we observe a close to exponential rise in the acceleration 
ratio as a function of global strain. The acceleration in- 
creases more slowly after strains of approximately five 
percent. 

This result also explains the observed changes to the 
persistence time distributions. In this case, the relevant 
parameter is the accumulated strain in between subse- 
quent hops. For instance, a strain step causes no change 
to p(r) since, after the first hop, there is no further strain. 
In the constant strain rate experiment, the acceleration 
factor can be defined in the same way as before, using in- 
stead the cumulatives of p(r). Plotted versus the strain, 
this gives the same universal curve as in Fig.^b) demon- 
strating that the second hops are accelerated in exactly 
the same way as the first hops. Although it is more diffi- 
cult to quantify, the dependence of p(r) on the wait time 
and the measurement time in the step stress experiment 
can be explained by the fact that the strain explicitly 
depends on these same variables. 



The universal collapse observed in Fig. |3jb) strongly 
suggests that the global strain is a good variable to de- 
scribe the influence of deformation on structural relax- 
ations, rather than the stress as postulated by the Eyring 
model. It is interesting to note that the strain was cho- 
sen as control variable in the Soft Glassy Rheology (SGR) 
model . Here the acceleration factor for barrier cross- 
ings increases proportional to exp(Z 2 ), where / measures 
the local strain on a mesoscopic element. In this model, 
the local strains increase in tandem with the global strain 
in between barrier crossings. Although the measured ac- 
celeration factor differs from that of the SGR model, in- 
spection of our particle trajectories shows that this may 
be a good picture of local strain in our system: individ- 
ual particles exhibit very little non-affine displacement 
except during a hop. 

We have presented the first investigation of the full 
relaxation time distribution of aging, deformed polymer 
glasses. As reported previously, all three deformation 
modes accelerate the segmental dynamics. Although the 
exact transformation of the relaxation time distributions 
is specific to the mode of deformation, we suggest that the 
acceleration can be described by a universal function of 
the global strain. Control simulations with longer chains 
containing 100 monomers indicate that this conclusion 
holds independent of chain length [TtJ • We also provide 
a possible explanation for the phenomenon of mechanical 
rejuvenation, whereby aging is arrested by deformation. 
The persistence time distribution causes aging in glasses 
because of the weak power law decay. Ergodicity is re- 
stored during flow under both a step stress and constant 
strain rate deformation because this distribution is nar- 
rowed and the mean persistence time becomes finite. 

This work was supported by the Natural Sciences and 
Engineering Research Council of Canada (NSERC). 



[1 
[2 
[3 
[4 

[5 

[6 

[7: 

[8 
[9 
[10 



Electronic address: jrottler@phas.ubc.ca 

L. S. Loo, R. E. Cohen, and K. K. Gleason, Science 288, 

116 (2000). 

A. F. Yee, R. J. Bankert, K. L. Ngai, and R. W. Rendell, 
J. Polym. Sci. Part Polym. Phys. 26, 2463 (1988). 
H.-N. Lee, K. Paeng, S. F. Swallen, and M. D. Ediger, 
J. Chem. Phys. 128, 134902 (2008). 

H.-N. Lee, K. Paeng, S. F. Swallen, and M. D. Ediger, 

Science 323, 231 (2009). 

H. Eyring, J. Chem. Phys. 4, 283 (1936). 

K. Chen and K. S. Schweizer, Europhys. Lett. 79, 26006 

(2007). 

A. V. Lyulin, B. Vorselaars, M. A. Mazo, N. K. Balabaev, 
and M. A. J. Michels, Europhys. Lett. 71, 618 (2005). 
F. M. Capaldi, M. C. Boyce, and G. C. Rutledge, Phys. 
Rev. Lett. 89, 175505 (2002). 

R. A. Riggleman, H.-N. Lee, M. D. Ediger, and J. J. 

de Pablo, Phys. Rev. Lett. 99, 215501 (2007). 

R. A. Riggleman, K. S. Schweizer, and J. J. de Pablo, 



■5 



Macromolecules 41, 4969 (2008). 
[11] L. C. E. Struik, Physical aging in amorphous polymers 

and other materials (Elsevier/North Holland, 1978). 
[12] A. V. Lyulin and M. A. J. Michcls, Phys. Rev. Lett. 99, 

085504 (2007). 

[13] H.-N. Lee, PhD thesis, University of Wisconsin (2009). 
[14] J. Baschnagel and F. Varnik, J. Phys.: Condens. Matter 

17, R851 (2005). 
[15] R. S. Hoy and M. O. Robbins, J. Polym. Sci. Part B: 

Polymer Physics 44, 3487 (2006). 



[16] M. Warren and J. Rottler, EPL 88, 58005 (2009). 

[17] M. Warren and J. Rottler, (to be published). 

[18] The acceleration factor can be accurately determined 
only when hop data is well resolved at times well be- 
low the cage escape time. For this reason, acceleration 
factors for only the longest t w are shown. 

[19] P. Sollich, F. Lequeux, P. Hebraud, and M. E. Cates, 
Phys. Rev. Lett. 78, 2020 (1997). 



